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Abstract. - We study the effect of spatial confinement on the strength of propulsive diffusio- 
phoretic forces acting on a particle that generates density gradients by exploiting the chemical 
free energy of its environment. Using a recently proposed simple model of a self-propelling device 
driven by chemical reactions taking place on some parts of its surface, we demonstrate that the 
force significantly increases in the presence of confining walls. We also show that such effects 
become even more pronounced in two-dimensional systems. 
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Last years have witnessed a growing technological, ex- 
perimental, and theoretical interest in scaling standard 
ihachinery down to micro- and nano-scales, such as pro- 
ducing pumps, motors, or sieves needed for the develop- 
ment of "lab on a chip" devices. For apphcations in, e.g., 
drug-delivery systems or micromechanics [1,2], one of the 
most challenging problems at this stage is to develop ways 
to enable small-scale objects to perform autonomous, con- 
trolled motion. Exploiting a variety of physical and chem- 
ical processes, a number of minimalistic, "proof of prin- 
ciple" proposals of such micro-engines have been put for- 
ward for both translational and rotational motion. Al- 
though the research in this area is still in its early stages, 
several of these proposals have already been tested exper- 
imentally (see, e.g., Refs. [3,4]; for a recent review, see, 
e.g., Ref. [2]). 

Among the key issues in the development of such en- 
gines are their "fuelling" and the tuning of the driving 
force and torque. Inspired by biological nano-engines that 
are using catalytic reactions to extract energy from the 
environment (e.g., kinesin moving along microtubules is 
a molecular motor that uses hydrolization of ATP as the 
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energy source), Whitesides and co-workers proposed a de- 
sign of self-propelling devices based on an asymmetric dec- 
oration of the surface of small objects by catalytic par- 
ticles promoting a chemical reaction in the surrounding 
medium [3]. The asymmetry in the catalyst placement 
leads to an asymmetry in the distribution of the reac- 
tion products, eventually also in the distribution of reac- 
tants, and this asymmetric distribution can provide motil- 
ity through a variety of mechanisms, such as surface ten- 
sion gradients, cyclic adsorption and desorption, and dif- 
fusiophoresis [2,5,6]; moreover, such a design removes the 
need of a "fuel reservoir" within the device because the 
reactants are provided by the environment. 

Following the idea of the asymmetric surface distri- 
bution of catalysts, Golestanian et al [7] have recently 
studied a model system in which the autonomous motion 
emerges as a result of self-created diffusiophoretic gradi- 
ents. Specifically, they considered a spherical particle of 
radius R with a single, point-like catalytic site fixed on its 
surface, which is immersed in the three-dimensional {3d) 
bulk of a reactive solvent. The catalytic site was deemed to 
promote within the solute a chemical reaction that creates 
product particles of a size much smaller than R. The spa- 
tial distribution of the emerging (diffusive) product parti- 
cles is asymmetric and thus results in an effective diffusio- 
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phoretic force acting on the particle. Several aspects of the 
resulting motion have been discussed, such as the ensuing 
transient time-dependent behavior, effects due to density 
and reaction-rate fluctuations, as well as the influence of 
rotational diffusion on the directionality [7] . A first exper- 
imental realization of this model, using platinum coated 
polystyrene spheres, has been recently reported. [8] 

For most of the applications in biological systems or 
'lab on a chip'-type devices one has to deal with a com- 
plicated internal structure of the system, e.g., networks of 
narrow channels or pores and various impenetrable impu- 
rities. Thus the spatial confinement is a relevant feature 
which may even lead to a quasi two- or one-dimensional 
behavior so that assuming the presence of an unconfined 
3d bulk reactive solvent, as in Refs. [1,2,7], may break 
down. Moreover, studies of bacteria or cell motion on 
surfaces such as nutrient substrates represent examples of 
two-dimensional systems. A stiidy of two-dimensional self- 
propelled particles with fluctuations in speed and direction 
of motion, aiming at describing such experimental systems 
has been recently reported. [9] Intuitively, one expects 
that the spatial confinement or the reduced dimension- 
ality may influence the resulting motion of self-propelling 
objects like the ones discussed above, but a priori it is not 
clear if this will lead to a decrease or an increase in the 
value of the diffusiophoretic force acting on the particle, 
and if such changes are of relevant magnitude. 

Based on the model used in Ref . [7] , here we address the 
effects of spatial confinement on the diffusiophoretic force 
exerted on a self-propelling particle. For a catalytic site 
with a uniform, i.e., time-independent reaction rate we 
show (1) that in 3d the presence of a confining wall for 
the reaction product leads to a significant increase of the 
force, even in the case of the simplest possible confining ge- 
ometry given by a concentric wall with the same shape as 
the particle, and (2) that a reduced dimensionality alone, 
i.e., a two-dimensional {2d) system, also strengthens the 
force acting on the particle. 

The system we consider is shown in Fig. 1. It consists of 
an impermeable, spherical particle (disk in 2d) of radius 
R with a point-like catalytic site (black dot in Fig. 1) on 
its surface (perimeter in 2d), which promotes a chemical 
conversion of a surrounding solvent into product particles 
(small hatched circles in Fig. 1). The particle and the 
surrounding solvent are enclosed in a concentric, imper- 
meable, spherical (circular in 2d) shell of radius Ri = rjR 
(r? > 1). Note that in 3d the model discussed in Ref. [7] 
is recovered in the limit 77 ^ oo. The whole system is 
assumed to be rigidly fixed at the common center O, such 
that we do not have to account for effects due to the mo- 
tion through the surroimding medium (e.g., if the solvent 
is a liquid, hydrodynamic interactions may also play a 
role), i.e., as a first step here we focus on the statics of the 
self-propelling device. 

Following Ref. [7], we assume that the reaction at the 
catalytic site acts effectively only as a point-like source of 
product molecules, which are diffusing in the solvent with 




Fig. 1: An impermeable, spherical (disk-like in 2d) particle of 
radius R with a point-like catalytic site at Vs (depicted as a 
black dot) on its surface (perimeter in 2d), enclosed by a con- 
centric, impermeable, spherical (circular in 2d) wall of radius 
Ri = rjR {rj > 1). The top inset shows the notations for the 2d 
case. The reaction products, indicated by small hatched circles, 
exhibit a gradient of their number density p{r,t) (per volume 
in 3d and per area in 2d). This setup can also be thought 
of as an elementary micropump creating a microcurrent into a 
confined environment. 



diffusion coefficient D; the amplitude of the production 
rate of this source, located at Vs = —Rz (where z is the 
unit vector of the z-axis), is denoted as A{t). Accord- 
ingly, the time evolution of the density p(r, t) of product 
molecules is governed by the diffusion equation 

dtp = DV^p + A{t)6^{r-r,), R < \r\ < Ri, (1) 

subject to the initial condition (IC) of zero density of prod- 
uct molecules and to the boundary conditions (BC) of zero 
normal current on the surface of the particle and on the 
confining wall: 

p(r,0) = 0, (2a) 

(^)l (-) 

As discussed in Refs. [7, 10], an asymmetric, non- 
uniform distribution of products p{r, t) induces particle 
motion with a velocity proportional to the first moment 
of the density distribution p{r,t) on the surface of the 
object. Since this is the result of a Stokes force-like cal- 
culation, we operationally define for our system, which is 
rigidly fixed, the net diffusiophoretic force acting on the 
particle at time t along the ^-direction as also being pro- 
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portional to the first moment of p{\r \ = R,t): 



while the regular part (i.e., the general solution of the 
homogeneous Helmholtz equation) can be written as 



F,{t) = -2nR^b3 j d0 sine cose p{R,0,t), in 3d, (3a) g{r,e,q;rs,v) = ^[a^(g; r^, r;)i^+i/2(gr) 
e>o 

Pe{q; Vs, r])ke+i/2{qr)]Pe{cose), 



(8) 



F^(t) = -Rh2 j cos (j)p{R,(l),t), in 2d, 



+ 



(3b) 



where 63 and 62 are some system-dependent constants with 
the dimension of energy and of the order of the thermal 
energy ksT [11] (in the following, we set 6 = 63 = 62 for 
simplicity). Knowledge of Fz{t) as a function of 77 < 00 
for confined 3d and 2d systems allows us to quantify the 
effect of confinement on the propulsive force, and thus we 
proceed by computing {t) for general values of r). 

The force is known once the density p{r,t) [i.e., the 
solution of Eq.(l) subject to the IC and BC conditions in 
Eq. (2)] is determined. Using the Laplace transform and 
the inversion theorem [12]: 
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(5) 



respectively, where f{p) is assumed to be well-defined 
for p G M+ and 7 G R is siifficiently large such that 
all singularities of / lie to the left of the integration 
path, one finds that the Laplace-transformed density is 
given by p{r,q) = —a{q)G{r,q), where q = \JpjD > 0, 
a(q) = A{(i)lD, and G(r,q) is the Green's bmction for 
the Helmholtz operator — with BCs of vanishing 
normal derivative at |r| = R,R\. Decomposing G{r,q) as 
G{r, q) = Gs{r, q) + g{r, q), where the singular part Gg is 
the free space Green's function for the Helmholtz opera- 
tor (which is known in any spatial dimension d, see, e.g., 
Ref. [13]), the initial problem has been reduced to that 
of finding the solution g of the homogeneous Helmholtz 
equation subject to the boundary conditions 



dg{r, q) 
dr 



\r\=R,Ri 



dGs 
dr 



(6) 
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While the main steps in the calculation of g are the 
same in 3d and in 2d, the details do depend on the 
dimensionality; for clarity, we will consider the two cases 
separately. 

The three-dimensional case 

In 3d the singular part of the Green's function is given 
by [13] 



Gs{r,q; r^) 



q e 



-q\r—rs \ 



An q\r — Vg 



(7) 



where^+i/2(z) = ^/TrJ{2z)Ii+i/2iz) and kf+i/2{z) = 
y^n/ {2z)Ki^i/2{z) are the modified spherical Bessel func- 
tions of the first and third kind, respectively [14], is the 
Legendre polynomial of degree £, while the functional co- 
efficients ai{q; rs,r]) and Peil', fs,v) will be fixed to fulfill 
the boundary conditions. 

Using one of the addition theorems for the Bessel func- 
tions [12], and noting that the angle between r and is 
TT — 9 (see Fig. 1), Gg can be re- written as 

Gs{r,q;rs) = - ^Q(g)i^+i/2(gr<)fc^+i/2(gr>)P^(cos6') 
e>o 

(9) 

where C£{q) — (—1)^(2/'+ l)q/{2n^), r< = min(r, r^), and 
r> = max(r, r^). By combining Eqs. (8, 9) and the bound- 
ary conditions [Eq. (6)], by noting that = i?+ (i.e., 
the source is on the spherical surface at r = i? so that 
— JS+ = 0+ ), and by re-writting the functional coef- 
ficients as ae = ceie+i/2{qR)c(e and f3e = Ciie+i/2{qR)Pe, 
the functions and /3£ are determined by the solution of 
the following system of linear equations: 



ae + Pewj,+i/2iqR) 
ai + /3ewi+i/2{r]qR) 



ve+i/2iqR), (10a) 

m+i/2{vQR), (lOb) 



where ^£+1/2(2) := A;£+i/2(z)/ii?+i/2(^^) and ^£+1/2(2) := 
[dk£j^i/2{z)/dz]/[di£^i^2{z)/dz]. This determines the 
Green's function G{r,q). 

Before proceeding, we make three remarks: 

(i) Since for qR > one has lim wgj^i/2{riqR) = 0, in 

?7 — *oo 

the limit r] ^ 00 one finds ag = and Eq. (10a) reduces 
to the corresponding BC in Ref. [7], i.e., as expected the 
solution for the unbounded case is recovered. 

(ii) Equations (10a) and (10b) will not coincide in the 
limit rj — > 1 because even in this limit a source singularity 
is present at and is picked up by Gs{q) [Eq.(9)]. 

(iii) The Laplace transform N{p) of the total number of 
product particles in the system at time t is 

R.l IT 

Nip = Dq^) = 2tt j drr^ j d6»sin6iPo(cos6»)[-a(j3)]G(g) 

R 

= -4na{p)coiq)ii/2{qR) 
Ri 

X J drr^[aoii/2iqr) + {Po-'i-)ki/2{qr)] (11) 

R 

t 

= Da{p)/p ^ N{t) = j dt'A{t') , 
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i.e., as expected N{t) is given by the time integral of the 
production rate, providing a welcome consistency check. 

We now focus on the particular case in which the ac- 
tivity of the catalytic site is time independent, i.e., A{t) 
is constant: A{t) = H{t)/Tf, where H{t) is the Heaviside 
step function and Tf is the average production time for 
the creation of a product molecule. Taking the Laplace 
transform of Eq. (3a) and noting that cos^ = Pi (cos 0), 
one obtains 



2'KhR^a{p) j d9sm9Pi{cos0)G{q) 





(12) 



where Fq := 2hR/{iTDTf) and 



+ [MVP\^)-'^]h/2{Vp)}- (13) 

At room temperature, using the order of magnitude es- 
timates from Ref. [7], i.e., 6 = fc^T, i? = 50 nm, 

= 5.5 X 10~^°m^/s (corresponding to a product-particle 
diameter ct = 0.8 nm and a viscosity of the solvent of 
H = I mPa X s), and 1/t/ = 25 kHz leads to values of 
Fq of the order of 10"^ pN, which would correspond to a 
velocity t;o ~ -Fb/ {l^R) of the order of 100 nm/s, while, for 
comparison, the weight W of such a particle made out of 
latex (density p ~ lO'^kg/m^, neglecting buoyancy) is of 
the order of 10^^ pN, i.e., Fq » W"- 

The complicated structure of the function $(p; rj) makes 
it rather laborious to carry out the inverse Laplace trans- 
form, so that the full time dependence of the force F^ 
cannot be derived easily. However, using the large p ex- 
pansion of Eq. (13) the short time behavior of F^ can be 
obtained from Eq. (12), 



F^{u < l,r])/Fo = VttV", u = Dt/R^ 



(14) 



This dependence can be rationalized as being due to the 
diffusion-like increase of the area on the particle around 
its catalytic site that is covered by product molecules. 

On the other hand, the asymptotic value Fz°°\'r]) := 
lim Fz{u;r]), i.e., the "steady state" force, can be deter- 

u—>oo 

mined using the inversion formula Eq. (5) by noticing that 

^{p; rj) has a simple pole at p = 0, which determines the 
asymptotic value as the residue of ^{p;r]) at p = [12]. 
This yields 



Fi°°)(7?)/Fo = -Res[$(p;r?),p = 0] 



which holds for any value of r; > 1. Together with Eqs. 
(12) and (15) the expansion of Eq. (13) for small p renders 
the asymptotic temporal approach of the force towards its 
steady-state value: 



F^(m> 1,7?) 
3 



t{7]) 



~ 1 — e""/'^'''^ , 



10 - 1 ■ 



(16) 



The corresponding relaxation time T{ri) = R'^/[DT{ri)] is 
finite even in the case of extreme confinement, i.e., for 
77 ^ 1: rij] ^ 1) = 1/2 and increases with decreasing 
confinement (increasing 77), as intuitively expected. 

For systems of infinite size (rj 00) we recover the 
resuh Fi°°^(r7 ^ oo)/Fo = 7r/4 in Ref. [7]. For all fi- 
nite T] one has Fz°°\ri) > Fz°°\f] — > 00). Therefore the 
spatial confinement leads to an enhancement of the dif- 
fusiophoretic propulsion, which is the main result of this 
analysis. As shown in Fig. 2, even at moderate values of 
rj the effect of the confinement is strong, e.g., at ~ 2 
there is a ca. 25% increase in the force compared with the 
unconfined case. For 77 — » 1 the force diverges; this is, of 
course, an artifact with stems from the assumption that 
the product molecules are point-like. Actually, there is a 
lower cut-off R1, where — i? is of the order of the hard 
core diameter a of the product molecules, below which 
this assumption breaks down and Eq. (15) is no longer 
valid. For cr <C i?, which is a reasonable assumption, one 
has 77c = l-|-o"/i?>l, and thus the steep increase in the 
force near 77 = 1 is physically relevant. 
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3(3?, ?/ ^ 00 
2d, ?/ ^ 00 




(15) 



7T 

1 

i] = R,/R 

Fig. 2: The steady-state force F^°°\rj)/Fo as a function of rj. 
The horizontal lines indicate the asymptotic values (7? — > 
oo)/Fq = 7r/4 in 3d and 7r/2 in 2d, respectively. 

The two-dimensional case 

The analysis follows along essentially the same lines as 
in 3d, so that we merely outline the steps involved and 
present the final result. Before proceding, we note that 
this case would correspond either to a very thin disk en- 
closed by a circular wall, or to a macroscopicaly long cylin- 
der with a longitudinal line of catalyst on its surface en- 
closed by another macroscopicaly long cylindrical wall. In 
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2d, the free space Green's function for the Helmholtz op- 
erator is given by [13] 

Gs{r,q;rs) = -^Ko{q\r-r,\), (17) 
and the regular part can be written as 

g{r,e,q;rs,v) = '^[am{q;rs,r])Im{qr) + (18) 

m>0 

+ l3rn{q;rs,'n)Km{qr)]cos{m(l)) . 

Here, Km{z) and Im{z) are the modified Bessel functions 
of integer index m [14], and the coefiicients am{q',rs,r]) 
and Pm{Q\fs,'n) will be fixed by the BCs. Note that in 
Eq. (18) there are no terms proportional to sin(m(/)) be- 
cause the solution has to be invariant with respect to the 
transformation (j)^ —<j). 

The addition theorem for the modified Bessel function 
Kq [12] together with the fact that the angle between r 
and Ts is equal to tt — </> (see the inset of Fig. 1) allows one 
to represent the singular part as the series 

{qr<)Km{qr>)cos{m(j)) (19) 

m>0 

where c„ = (-1)™(2 - (5„,o)/(27r). Using Eqs. (18), (19), 
the BCs [Eqs. (6)], and noting that rg = R+, the coeffi- 
cients Orn = CrnIrn{qR)am and = CrnLm{qR)Prn follow 

from the solution of the system of linear equations 



WmiqR) 



v-miqR) , 

WrnimR) : 



(20a) 
(20b) 



where Vmiz) := Km{z)/Im{z) and Wm{z) := 
[dKm{z)/dz]/[dIm{z)/dz]. This fixes the Green's 
function G{r, q). 

Similarly to the 3d case, in the limit r/ ^ oo one recovers 
the result of a direct calculation for the unconfined system, 
the limit rj ^ I still picks up the singularity at r,,, and 
the total number of particles is given by the time integral 
of the production rate. 

Again focusing on the case of a time- independent source, 
A{t) = H{t)/Tf, from the Laplace-transformed Eq. (3b) 
one obtains 

27r 



F,ip) = bRaip) J dcP{cos<l>)G{q) 



(21) 



Dt 

i?2 



where 



-|(£-M$i(p;r/)]) \u, 



$i(p;ry) = ^-^-^^ x [ai{y/p;ri)h{y/p) 

+ ms/p;r,)-l)Ki{^)\ . (22) 



The short time behavior of Fz is obtained from Eq. 
(21) by using the large p expansion of Eq. (22). One 
finds the same expression [Eq. (14)] as in 3d, which sim- 
ply refiects the increase over time in the "support" of the 
distribution of the product molecules on the surface of 
the particle. The "steady state" diffusiophoretic force 
Fz°°\r]) := lim Fz{u;rj) follows from the inversion for- 

u — *oo _ 

mula in Eq. (5) and the fact that ^i{p;r]) has a simple 
pole at p = 0: 



^ z 

Fi~)(r?)/Fo 



{n)/Fo = --Res[$(p;77),p=0] 
nr]^ + 1 



2 772 



1 



(23) 



The long time behavior, i.e., the approach of F^ to its 
asymptotic value follows from Eqs. (21) and (23) by using 
the small p expansion in Eq. (22) leading to 



Fz{u > l,r?) 

iv) 



p(oo) 



1 



-u/T{ri) 



1 rj^ ln(77) 



1 



(24) 



F^'^'iri 



The relaxation time T{ri) = R^ /[DT{ri)] stays finite even 
for extreme confinement because t(?7 ^ 1) = 1/4, and it 
increases with decreasing confinement, but only as ln(77), 
i.e., much slower than the increase ^ 77^ in the 3d case (see 
Eq.(16)). We note that occurrence of terms ~ sln(s) or 
~ t/ln{t) in the Laplace or the time domain, respectively, 
is generic for two-dimensional infinitely large reaction- 
diffusion systems [15] because the fractal dimension of 
Brownian paths coincides with the dimension of the em- 
bedding space. For confined systems, however, the loga- 
rithmic correction ln(t) to the relaxation time is replaced 
by ln(?7), which is related to finite-size effects. 

Equation (23) implies that the value of the force 

00) /Fq = 7r/2 in the unbounded case is twice as 
large as the corresponding one in 3d (assuming the same 
values for R, D, b, and r/ in both cases). For any rj, 
the force Fz°°\'r]) is stronger in 2d than in 3d, which can 
be easily inferred from the ratio of Eqs. (23) and (15), 
as shown in Fig. 2. Thus we conclude that reduced di- 
mensionality enhances the propulsion mechanisms. Such 
a two-dimensional system can be realized, e.g., by using a 
disk-like particle floating on a liquid-vapor interface and 
product particles that are mobile on the interface but do 
not dissolve in the bulk liquid and do not evaporate either. 

As in the 3d case, with increasing confinement (i.e., 
decreasing 77) the force increases steeply from the value 
Fz°°\r] 00)/ Fq = Tr/2 and for 77 ^ 1 it diverges (up to 
a lower cut-off at ?7c > 1 + a/R). 

In summary, based on the model proposed in Ref. [7] 
we have studied the effect of confinement on the diffusio- 
phoretic force exerted on a self-propelling particle in both 
three and two spatial dimensions. The main findings (see 
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Eqs. (15), (23) and Fig. 2) are: (1) The presence of a con- 
fining wall for the reaction product leads to a significant 
increase in the force acting on the particle, even in the case 
of the simple concentric confining geometry (see Fig. 1). 
(2) Reduced dimensionality alone (i.e., even without con- 
finement) is already sufficient to increase the efficiency of 
the proposed propulsion mechanism, leading to a stronger 
force on the particle. A possible way for testing such pre- 
dictions would be to use a colloidal particle with a catalyst 
grafted on its surface and enclosed in shape-controllable 
vesicles filled with a solution of reactants. By superposi- 
tion the approach presented here can be extended to the 
case of an arbitrary decoration of the particle surface by 
catalytic sites as long as the reactions can be assumed to 
occur independently, i.e., without competing for the reac- 
tants in the common environment. This will allow one to 
find an optimal decoration for providing stability against 
particle rotations, which can otherwise spoil the imidi- 
rectional motion. Further extensions may focus on new 
phenomena emerging from more complicated geometries. 

* * * 

M.N. P. gratefully acknowledges very fruitful discussions 
with A. Gambassi, L. Harnau, and M. Tasinkevych. 
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